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The First Reported Infrared Emission from the SN 1006 Remnant 

P. Frank Winkler^, Brian J. Williams^, William P. Blair^ Kazimierz J. Borkowski^, Parviz 
Ghavamian^, Knox S. Long^, John C. Raymond^, and Stephen P. Reynolds^ 

ABSTRACT 

We report results of infrared imaging and spectroscopic observations of the 
SN 1006 remnant, carried out with the Spitzer Space Telescope. The 24 /im image 
from MIPS clearly shows faint filamentary emission along the northwest rim of 
the remnant shell, nearly coincident with the Balmer filaments that delineate the 
present position of the expanding shock. The 24 /xm emission traces the Balmer 
filaments almost perfectly, but lies a few arcsec within, indicating an origin in 
interstellar dust heated by the shock. Subsequent decline in the IR behind the 
shock is presumably due largely to grain destruction through sputtering. The 
emission drops far more rapidly than current models predict, however, even for 
a higher proportion of small grains than would be found closer to the Galactic 
plane. The rapid drop may result in part from a grain density that has always 
been lower — a relic effect from an earlier epoch when the shock was encountering 
a lower density — but higher grain destruction rates still seem to be required. 
Spectra from three positions along the NW filament from the IRS instrument 
all show only a featureless continuum, consistent with thermal emission from 
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warm dust. The dust-to-gas mass ratio in the pre-shock interstellar medium is 
lower than that expected for the Galactic ISM — as has also been observed in the 
analysis of IR emission from other SNRs but whose cause remains unclear. As 
with other SN la remnants, SN 1006 shows no evidence for dust grain formation 
in the supernova ejecta. 

Subject headings: ISM: individual (SNR SN1006) — ISM: kinematics and dy- 
namics — shock waves — supernova remnants 



1. Introduction 



The supernova of 1006 CE. is generally recognized as the brightest stellar event in 
recorded human history, with surviving contemporary reports from China, Japan, Korea, 
the Arab world, and in Europe as far north as St. Gallen, Switzerland (latitude 47 . 4° N) , 
all despite the supernova's southern declination of Decl. (1006) = —38.5° (IStephensonI |20 1 Ol ) . 
Throughout much of the twentieth century, SN 1006 languished in relative obscurity com- 
pared with the attention lavished on th e other eleventh -century supernova, SN 1054, that 
gave rise to the Crab Nebula. But once iGoldsteinl ( 119651 ) called attention to the event and 
MarsdenI (119651) constrained its l ocatio n, radio emission from the SN 1006 remnant was soon 
discovered by [Gardner fc Milnd (Il965[ ). It has si nce become a wid e ly obs erved and important 
object: X-rays f rom S N 1006 were discovered by I Winkler &: LairdI ( 119761 ): optical emission by 
van den Berghl ( 1l976l ) ; and (following a somewhat checker ed history) TeV 7- r ays were defini- 
tively detected and mapped by the HESS collaboration (lAcero et al.l 120101 ). iKoyama et al. 
(119951 ) discovered that the NE and SW limbs of SN 1006 have a featureless, power-law X-ray 
spectrum, providing the first conclusive evidence that electrons can be accelerated to TeV en- 
ergies in supernova shocks and cementing the long-suspected connection between supernovae 
and cosmic rays. Yet no infrared emission has previously been reported from SN 1006. 

The SN 1006 remnant has a limb-brightened shell, 30' in diameter with st rong bi-lateral 
symm etry in radio and X-rays. Its distance is well determined at 2.2 kpc (jWinkler et al. 
20031 ). and it is located high above the Galactic plane {bu = 14.6°, z = 550 pc). Its high- 
latitude location leads to relatively low foreground absorption. SN 1006 itself was almost 
certainly a Type la event, based on its location, environment far from any recent star forma- 
tion, apparent absence of any compa ct remnant, and implicatio n from Chinese records that 
it remained visible for several years ( IStephenson fc Clarkll2002l ). The optical emission frorn 



its remnant consists solely of Balmer-domin ated filaments (e.g., iSchweizer fc Laskerl Il978 
Kirshner et al.l Il987l : iGhavamian et al.l |2002|) , which indicate that it is expanding into an 
environment with relatively low density that is at least partially neutral, another signature 
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of SN la remnants. 



Supernovae and their remnants are critical to dust grain formation, modification, and 
destruction. Dust grains may condense from metal-ricli SN ejecta, and the fast shock waves 
in SNRs can lead to sputtering thr ough collisions with shock-heated gas that erode an d 
eventually destroy dust grains (e.g., iDwek fc Scald Il980t iDraind l2003l : ISankrit et al.l 120101 ). 
Although these processes have been observed in a number of SNRs, SN 1006 nevertheless 
provides a unique perspective due to its location high above the Galactic plane in a relatively 
pristine region of the ISM. There is growing observationa l evidence for the formation of dust 



in th e ejecta from core-collap se SNe, i ncluding SN1987A (iWooden et al.lll993l : lde Kool et al. 



Nebula (I Gomez et al 



2012b 



SNe, (e.g.. iKotak et al 



19981 : iMatsuura et al.ll201lh . Gas A fiRho et all l2008l : iNozawa et al.l 12010^ . and the Grab 



2009 



, and in the late-time emission from a number of more distant 



Szalai et al.ll201ll : iMeikle et al.ll201l[ ) 



Recent models by iNozawa et al.l (120111 ) predict that small dust grains can also form in 
the ejecta from Type la SNe, and that 1O~^M0 of dust can survive for ~ 1000 years. How- 
ever, no dust associated with ejecta has as yet been observed in any of historical remnants 
of the SNIa: Tycho (SN1572), Kepler SN1604), or RGW86 (probable SN185 remnant), 
all of which have been widely observed in the infrared. Mid-IR emission f rom both Ty - 



20071 : IWilhams et al.ll2012l ). Herschd 



was also detected from IRAS f lArendtl Il989l : ISaken et al 



Gomez et al. 



20 



rom ISO (iDouvion et al.l 1200 ll ). Spitzer (iBlair et al. 



cho and Kepler was detected from the InfraRed Astronomy Satellite (IRAS, iBraun 
and both have been observed in detail 



1987) 



2a), and other missions. RGW86 



1992 



andlWilliams et al.l teOllaf) 



Borkowski et al. ( 20061 ) have 



have recently presented detailed observations with Spitzer. 
also reported IR observations for four SNIa remnants in the Large Magellanic Gloud. All 
these cases show thermal emission from dust grains and provide evidence for grain destruc- 
tion in SN shocks, but none shows convincing evidence for grain formation (see especially 
Gomez et aDl2012ah . 



In this paper we report the first detection of IR emission from the SN 1006 remnant. In 
Sections 2 and 3 we describe the results of an imaging mapping campaign with the Spitzer 
Space Telescope, which reveals a filament of 24 fim emission along the NW limb of the 
shell, plus very faint, diffuse emission elsewhere. In Section 4 we describe the spectroscopic 
observations from the InfraRed Spectrometer (IRS) on Spitzer, and in Section 5 we present 
a discussion in which we fit dust models to both sets of data and assess the total dust mass 
and dust-to-gas mass ratio. Finally, we give a brief summary of our results in Section 6. 
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2. MIPS Observations 

We mapped SN 1006 twice witli tlie Multiband Imaging Pliotometer for Spitzer (MIPS, 



Rieke et al.ll2004( ) as part of program 30673: on 2007 March 5 and 2007 August 22. In order 
to cover the entire 30' shell and the surrounding background field, each observation required a 
pair of overlapping scan maps, which were taken consecutively. The March AORs (18724608 
and 18725120) were carried out during a "cold" campaign and included observations at 
24/im, 70//m, and 160/im, while the August ones (18724864 and 18725376) were during a 
"warm" phase of the cryogenic mission and gave only 24/im and 70/im data. Both pairs 
of scans consisted of 14 legs of length 0.75°; offsets of 148" were used for the March pair, 
while the August pair used offsets of 111" to give greater exposure of SN 1006 itself at a cost 
of somewhat less surrounding background. We began with Basic Calibrated Data (BCD) 
files from the Spitzer Heritage Archive that had been processed through version 18.12 of 
the Spitzer pipeline, and we then performed further fiat-fielding by taking a median of the 
hundreds of frames from each AOR, and dividing all the frames from that AOR by the 
normalized flat. The estimated zodiacal light was then subtracted from each frame, and the 
maps were mosaicked together onto a standard world coordinate system, 50' square with 1" 
pixels, using the MOPEX software provided by the Spitzer Science Center. Finally, we used 
DAOPHOT in irafQ to remove the vast majority of point sources (mainly stars) from the 
mosaic image, in order to better assess the diffuse emission from SN 1006 itself. 

The full 24yum mosaic is shown in Fig. 1, with both the original and star-subtracted 
versions in panels a and b, respectively. These clearly show the first detection of SN 1006 
in the mid-infrared. For reference, we also show images of the identical SN 1006 field in 
Ha, X-rays, and radio in other panels of the figure. The most prominent IR feature is 
the thin filament along the NW rim of the shell. Even this brightest emission is quite 
faint, with peak surface brightness levels above the background of 0.25-0.3 MJy sr~^, only 
about 6% of the background emission and fainter than either the Tycho or Kepler SNRs 
by more than a factor of 20. The integrated flux from this filament in the 24/im band is 
F(24/im) = 130 ± 30 mJy, where the precision is limited by our tracing of the filament and 
ability to assess the background accurately, since there is diffuse background as well as faint 
residual stars in the "star-subtracted" mosaic. 

We used a similar approach for the 70 fim data, except that here the mosaic has a 
scale of 4"pixel~^. Despite our flat-fielding all the data using median frames for each AOR, 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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artifacts along the scan directions are nevertheless evident in the final mosaic, shown in the 
bottom left panel of Fig. 1. There is the slightest hint of an enhancement in the region 
of the NW filament, but this does not trace the narrow filament and is barely larger than 
background fiuctuations in the vicinity. We give only an upper limit for the average surface 
brightness in the region of the filament: S'(70/im) < 0.85 MJy sr~^ (3ct), and an integrated 
fiux F(70/im) < 800 mJy. The 160 fj,m images, compiled using only the two scans from 2007 
March for which 160 fxm are available, are quite noisy and show no discernible emission from 
SN1006. 



3. Imaging Analysis and Multi-wavelength Comparisons 

3.1. The Northwest Filament 

The NW filament c orresponds closely with the optical filament that appears only in the 
hydrogen B aimer lines (Ivan den Berghlll976l : IWinkler et al.ll2003l : [Raymond et al.l 120071 ). as 
shown in the deep continuum-subtracted image shown in the lower right panel of Fig. 1. 
Such "nonradiative" filaments characterize fast shocks expanding into a low-density, par- 
tially neutral environment, and that are generally associated with relatively young remnants 
from Type la supernovae such as SN 1006. Balmer emission from these filaments has two 
components, both the result of neutral H atoms that cruise through the shock. In the hot 
post-shock environment, the neutrals can be excited and decay to produce narrow lines (with 
a width characteristic of the pre-shock temperature), or they can undergo charge exchange 
with hot protons to prod uce broad lines whose width is closely related to the pos t-shock 
proton temperature (e.g., IChevalier et al.lll980l : iGhavamian et al.l 120021 : lHeng||2010l ). Since 
the lifetime of neutral atoms in the hot post-shock environment is very short, the Balmer fil- 
aments can occur only immediately behind the shock, and thus delineate the current position 
of the shock. 



To investigate the post-shock characteristics of emission along the NW filament in vari- 
ous bands, we have compared 24/im, Ha, and X-ray images. We have used a higher resolution 
Ha image, taken from the 0.9 m telescope at CTIO on 2002 March 2 1, for which the star s 
were subtracted using a matched continuum image (for details, see IWinkler et al.l l2003l ). 
instead of the one shown in Fig. Ic. Likewise, in X-rays we have used the 0.5-2.0 keV data 
from a 90 ks obs ervation taken 2001 April 26-27 with the Chandra ACIS-S (for details, see 
Long et al.l 120031 ). rather than the mosaic of shorter ACIS-I exposures used in Fig. le. We 
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first checked the ahgnment of all three imaged using four near point sources common to all: 
one bright star and three distant galaxies, and found their positions in all three bands to 
agree within 0.5". To simulate images in all three bands at a common epoch, we used the 
/i = 280 ± 8 mas yr~^ proper motion rate for the NW Balmer filaments, measured relative 
to the center of cur vature at R.A. (2000.) = 15^ 03°^ 16!2, DecL (2000.) = -42° 00' 30" by 



Winkler et al.l (120031 ). to extrapolate the Ha and X-ray images from their original epochs 
(2002.2 and 2001.3, respectively) to the 2007.5 mean epoch for the SpitzerMlPS scans. (Re- 
cent measurements from Chandra find X-ray proper motions alon g most of the NW r im of 



SN 1006 that agree with those measured for the optical filaments, iKatsuda et al.ll2012l ). 



In Fig. |2]we show a three-color combination that compares the 24 /xm (red). Ha (green), 
and X-ray (blue) images of the NW filament, and in Fig. E] we plot radial profiles taken 
perpendicular to the filament, integrated over the four azimuthal regions outlined in Fig. 
[21 It is clear that the 24 /xm emission has a broader spatial profile that peaks 3-5 arcsec 
behind the Balmer filament, as best illustrated in region D of Fig. [3]where the profiles are the 
cleanest. The filament is also noticeably broader in the IR than in Ha, with FWHM > 15". 
This is not simply a resolution effect. The resolution of the 24 fim mosaic image — measured 
from profiles of several isolated stars — is 6" (FWHM), consistent with the diffraction limit, 
but small compared with the width of the filament. 

We attribute this IR emission to dust from the ISM that enters the shock, where it 
is heated through collisions with hot post-shock gas, produces IR emission, and is then 
destroyed. Somewhat further behind the shock, soft, primarily thermal. X-ray emission 
appears, rising to its peak about 10" behind the shock, just as the 24 /im emission is decaying 
away. This indicates the longer time scale for pro ducing the highl y ionized species that 



radiate efficiently and dominate the X-ray emission (ILong et al.l 120031 ). We can estimate the 
time scales through simple kinematics: the primary shock is moving outward to the NW at 
fi = 280 ± 8 masyr"^, and if we assume the nominal compression ratio of 4 between the 
post-shock and pre-shock densities, the post-shock material is dropping behind the shock at 
speed fs/4, or /i = 70 mas yr~^ (relative to the shock). Therefore, the IR peak occurs some 
60 yr after that material was shocked, and the X-ray peak occurs ~ 150 yr after the shock. 
We shall discuss the emission profiles further in Sec. 5.1 of this paper. 

In 2012 March , all-sky survey data from the Wide-field Infrared Survey Explorer (WISE, 
Wright et al. 2010) were released|f| Examination of data from the SN1006 region in all the 



^For checking the Ha image we of course used the version before continuum subtraction; the coordinates 
for the subtracted image are identical. 

see |http : / / wise2 . ipac.caltech.edu /docs/ release / allsky / 
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WISE bands (3.4, 4.6, 12 and 22 fim) shows emission from SN 1006 only in the 22 fim 
band. Not surprisingly, the WISE 22 /im image is very similar to the Spitzer 24 fim one, 
and clearly shows the NW filament. However, the WISE data have somewhat lower angular 
resolution and signal-to-noise than that from Spitzer. The non-detection of emission in the 
shorter-wavelength WISE b ands is consistent wit h the lack of short-wavelength emission 
in other Type la remnants (IBorkowski et al.ll2006l ). Cl ose examination of the images frorn 
the Improved Reprocessing of the IRAS Survey (IRIS; iMiville-Deschenes fc Lagachd l2005l ) 
also shows traces of the NW filament in the co-added 25 /^m image; however, it is unlikely 
that this faint feature would ever have been independently discovered in these low-resolution 
(1.5') maps. 



3.2. Other Faint IR Emission 



In addition to the narrow filaments along the NW rim of the shell, there is patchy, 
diffuse emission within much of the southern part of the SN 1006 shell, clearly visible in the 
24 /im mosaic and less so in the 70 fim one. The morphology of this emission is similar at 
both 24 /im and 70 fim, and is also similar to patchy emission extending well outside the 
SN1006 shell, especially to the south and east. In the southern portion of SN1006, there 
is also faint, diffuse emission in both Ha and X-rays (Fig. Ic and e), but unlike along the 
NW rim, the morphology in the south is only vaguely similar to that seen in the IR. In the 
absence of further data, we cannot determine whether the faint IR emission is associated 
with SN1006 itself or stems from cool foreground dust. 

In X-rays, SN 1006 shows knots of emission interior to the shock that are thought to 
arise largely from SN ejecta, but there are no obvious small-scale IR features that correlate 
with these X-ray-emitting ejecta knots. Hence we find no evidence for the formation of dust 
in SN 1006 ejecta. This is somewhat surprisin g in view of theoretical predictions for dust 
formation in S NIa electa (jNozawa et al.ll201lf), but is cons istent with observations of other 



Type la SNRs (IBorkowski et al. 



2006 



WiUiams et al.ll2012l ) 



4. Spectroscopic Observations 

We conducted IRS observations, all in Staring mode for about 90 min each, at three 
locations along the NW filament. We show the slit positions of these IRS pointings, labeled 
1-3, in Fig. SI overlaid on the MIPS 24 fim image of the filament. As noted earlier, the 
filament is quite faint in the IR (its average surface brightness is only about 6% above the 
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background at 24 /iin), and it was undetected in the SL (5.5-14 fim) spectrum. However, 
faint emission is detected in the LL order 2 (14-19.3 fim) spectrum, where the dust continuum 
begins around 17 /im and continues to rise throughout the LL order 1 (19-40 /xm) spectrum. 
We see no evidence for hne emission anywhere in the spec trum, consis t ent with the non- 



radiative nature of th e blast wave shock at this location (IDwekl 119871 : iDwek et al.l Il996 
Wilhams et aPboilbh . 



The data used in our analysis (AOR IDs 28138240 and 24142080) were pr ocessed with 



Spitz er calibration pipeline version S18.7.0. Using the IDL CUBISM software (jSmith et al. 



20071 ). we assembled the 72 BCD files from positions 1 and 2 along the NW filament into two 
data cubes, one for LLl and one for LL2. The IRS spectrum from position 3 is similar, but 
we have omitted it from our subsequent analysis becau se it falls on a re gion where different 
abundances were required to fit the X-ray spectrum (ILong et al.ll2003l ). Hot pixels affect 
the BCDs (especially beyond 35 /im), so we median smoothed the BCDs before creating the 
data cubes using the FILTERJMAGE routine from the IDL Astronomy User's library. The 
smoothing replaces the value of each pixel with the median of the surrounding 3 pixels in 
a moving box and significantly mitigates the hot pixels in the resulting data cubes. While 
the smoothing reduces the spatial resolution of the data somewhat, the emission from the 
filament is highly uniform and we expect minimal loss of information. 

After combining the BCDs within CUBISM, we extracted LLl and LL2 spectra from 
the NW filament using positions 1 and 2 marked in Fig HI These two positions are only 
~ 30" apart and their IRS spectra appear virtually identical; therefore, we have combined 
them into a single object spectrum. To estimate the sky background, we extracted separate 
LLl and LL2 spectra from locations just outside SN 1006 (marked in FigH]). We averaged 
the two sky spectra in each order and combined them to produce a final sky spectrum. We 
subtracted this spectrum from that of the SN 1006 filament to obtain the final spectrum 
shown in Fig. [5l We describe the fits to this spectrum in the following section. 



5. Discussion 

5.1. Model Fits to the Data 

Infrared emission in typical SNR shocks is produced entirely by collisional heating of 
dust grains (Dwek & Arendt 1992), where energetic electrons and ions collide with the grains, 
heating them to temperatures of 50-100 K. Collisions with protons and heavier ions are also 
responsible for grain sputtering, reducing the size of large grains and destroying many small 
grains entirely. We have developed models for such dust that self-consistently calculate both 
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the heating and the destruction of grains in the post-shock gas as a function of the density 
of the hot p lasma, rip, the electro n and ion temperatures, and the shock sputtering age 
r = J^Updi ( Borkowski et al. 2006 ). For an arbitrary input grain size distribution, we can 
use these models to calculate the IR spectrum and the post-shock profile that results as 
grains are heated, and eventually destroyed. 

The electron temperature, Tg, and the ionization tir nescale can be det ermined from a 
non-equilibrium ionization analysis of the X-ray spectra, as lLong et al.l (120031 ) did for the NW 
limb using Chandra data. We have revisited that analysis, taking into account hydrocarbon 
contamination of the ACIS detectors that was unrecognized at the time of the original 
analysis, and using more reliable atomic data. We first redetermined the absorption column 
density N^ from the reprocessed sy nchrotron- dora inated eastern limb data (region NE-1 of 
Long et al.ll2003[ ). fit using XSPEC (jArnaudlll996l ) with a pure synchrotron srcut model to 



obtain Nb = (8.2 ± 0.8) x 10 cm~ . We fixed A^h at this value in fitting Long et al's region 
NW-1, which includes our IRS regions 1 and 2. We obtained reasonable fits using a solar- 
abundance plane-shock model (vnpshock) , and obtained values of kTe = 0.9 keV, ionization 
timescale n^t = 5.3 x 10^ cm~^ s = 170 cm~^ yr, and emission measure n^Mg = 6.2 x 10~^ Mq 
cm~^. This is a s omewhat higher electron temperature than the kT^ = 0.6 keV found by 
Long et al.l (|2003[ ). but is still far below the post-shock proton tempe ra ture of kTp = 16 



keV, inferred from the Balmer line profiles by iGhavamian et al.l ( 120021 ). (IVink et al.l 12003 



also measured a similar value for Tg at a different position along the NW rim.). We then 
applied a Coulomb- heating model for collisions between electrons and protons behind the 
shock, and found that over the relatively small region we consider here, the mean proton 
and electron temperatures are 16 keV and 1.0 keV, respectively. We use these temperatures 
for the remainder of this work. 

Our models also require as input a distribution of grain sizes in the ISM ahead of 
the sh ock. A commonly used distribution for the Galaxy is that of IWeingartner fc Draine 
( 120011 ). Their models are appropriate for the diffuse ISM in the disk of the Galaxy, where 
the total-to-selective extinction ratio Ry = Av/E{B — V) = 3.1, and for dense clouds where 
Rv > 3.1. At the location of SN1006 some 500 pc above the Galactic plane, however, 
the abundance of dust — both in total amount and in the distribution of grain sizes — is likely 
very different from that near the Galacti c plane. There is signific a nt evidence for much lower 
values for Ry at high Galactic latitudes (ILarson fc Whittetll2005l ). iMazzei fc Barbarol (l201ll) 
have c alculated grain size distributions using the same formalism as IWeingartner fc Draine 
( 1200 ll ) for over 70 lines of sight to stars with "anomalous extinction sightlines." Particularly 
relevant for our work, they extend these mode ls down to low-density e nvironments with 
Rv = 2.0. We use the model from Table 3 of iMazzei fc Barbarol koil\ \. with Ry = 2.0 
and 10^6c = 5.0 (this is a parameterization in their models that corresponds to the number 
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of small carbonaceous grains). The grain size distributions corresponding to lower values 
of Ry are steeper than their "typical" Milky Way counterparts, i.e., they have an excess of 
small grains relative to large ones, resulting in a larger fraction of the total mass in grains 
being destroyed via sputtering for a given shock, as well as a different gas density required 
to produce the observed spectrum. 

In our models, the spectrum from warm dust is most sensitive to the gas density in 
post-shock plasma, so we have used proton density Up as the principal free parameter. Using 
the grain size distribution described above, we have fit the IRS spectrum over the range 21 - 
37 fim and obtain a best fit for a density of = 1.4 ±0.5 cm~^, shown in blue in Fig. [51 The 
21 - 37 fim range comprises all the data for which the signal-to-noise ratio is high enough 
to be useful. The uncertainty quoted is a formal 90%-confidence limit based on statistical 
uncertainties only, determined from the rms deviations from an arbitrary smooth polynomial 
fit to the spectrum in this range. We have made no provision for systematic errors in the 
IRS spectra, which are difficult to quantify. 

There is good reason to believe that the actual density is near the lower limit of the 
fitted range, near Up = 1 cm~^. From an analysis of the th ermal X-ray emission along 



the NW rim based on Chandra observations, iLong et al.l (120031 ) found a post-shock electron 
density Ue ~ 1 cm~^. Ana l ysis o f X-ray data from XMM-Newton using a different modeling 



approach led lAcero et al.l (120071 ) to a pre-shock density value of uq ~ 0.15 cm~^, equivalent 
to a post-shock density of Up ^ 0.6 cm~^ for a standard compression ration of 4. And finally, 
an estimate based on the thi c kness of the Ha filaments in an imag e from the Hubble Space 



Telescope by [Raymond et al.l ( 120071 ) , and subsequently corrected by iHeng et al.l ( 120071 ) , gave 
no = 0.15 — 0.3, and thus a post-shock density np ~ 1 cm~^. In Fig. |5] we also show, in red, 
a model with Hp = 1.0 cm~^, which gives an acceptable fit to the IRS data. Finally, we note 
that the predicted 70 fim surface brightness for any of our models is < 0.25 MJy sr~^, well 
below the upper limit reported in Section 2. 

We have also calculated the spatial profile for the decay of IR emission as grains are 
destroyed through sputtering in the hot, post-shock environment, assuming that the shock 
has been advancing through an ambient medium of constant density. In comparing the 
model with the observed radial profiles, we find that the observed decay is far more rapid 
than the sputtering model predicts, as illustrated in Fig. |6l Even if we assume a grain- 
size distribution heavily weighted toward small grains to give Ry = 2.0, for a density of 
nil = 1 cm~^ the sputtering rate would need to be artificially enhanced by a factor 6 to 8 to 
match the observed profiles. 

However, it is almost certain that the shock along the NW rim of SN 1006 has not been 
passing through a uniform-density environment. The ambient density to the NW is clearly 
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higher than that around the remainder of the shell, as evidenced by the following facts: (1) 
the flattening of the shell to the NW; (2) the fact that the proper motion a long most of the 
NW limb is only about 60% as high as that to the NE ( Katsuda et al. 12012); (3 ) the density 
contr ast around the shell required from X-ray observations ( Long et al.l 2003 : Acero et al. 
20071 ): and (4) the mere fact that there is significant IR emission along the NW limb, yet 
little or none elsewhere. Since the shell radius to the NW is only about 15% smaller than 
elsewhere, the shock must have encountered the denser region there relatively recently. As 
a result, the post-shock density profile, for both gas and dust, should drop as one moves 
inward from the shock. Therefore the decrease in IR emission behind the shock is likely 
due in large part to the fact that farther inward the grain density has always been lower, 
representing an earlier epoch in the evolution of SN 1006. We suggest that this may account 
for the rapid drop in IR emission ~ 20" — 35" behind the shock (Fig. For a density 
Up = 1 cm~'^ (equivalent to rig = 1.2 cm~^), the shock age from X-ray spectroscopy is about 
170/ne = 140 years, suggesting that the shock was encountering considerably lower densities 
earlier. Detailed modeling of such a scenario is beyond the scope of the present paper. 

Nevertheless, the fact that the observed post-shock drop in IR emission is so much 
sharper than predicted by the models, even close to the shock, probably indicates that 
current models significantly underestimate the actual dust sputtering rate. Effects that could 
be responsible include an even greater proportion of small grains, grains of high porosity, or 
a grain composition that leads to higher sputt e ring r ates. Enhanced sputtering rates are not 
without precedent; ISerra Diaz-Cano &: JonesI ( 120081 ) found that hydrogenated dust grains 
(perhaps the most common type of grains in the ISM, at least for carbonaceous materials) 
are sputtered at rates several times higher than their pure amorphous counterparts. The 
spectral data are insufficient to determine whether the carbonaceous grains in SN 1006 are 
hydrogenated or amorphous. 



5.2. Dust Mass of the NW Filament 



Almost independent of the models, we can determine the total amount of dust radiating 
in the NW filament, provided that the integrated IR luminosity is relatively well-known. We 
have used the measured IR fiux, F(24/im) = 130 ± 30 mJy (Section 2), and assumed that 
the overall shape of the dust spectrum observed with IRS (shown in Fig. E]) is applicable 
along the entire length of the filament to obtain Lm = (1.7 ± 0.9) x 10^^ ergs s~^, where 
we have scaled to a distance of 2. 2 kpc: ^2.2 = d/2.2 kpc. Our models assume optical 
constants from iDraine &: Led ( 119841 ) for graphite and amorp hous silicate grains i n the ISM, 
with relative proportions and grain-size distributions from iMazzei &: Barbard (1201 if ). As 
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m 



Williams et al.l (|2008[ ). we then calculate the dust mass in the NW filament as Mjust 



(1.1 ± 0.6) X 10~^(i2^2 where the uncertainty limits for both Lir and Mjust include the 
overall 15% uncertainty in the SpitzerMlVS extended source flux calibration. Of course, this 
method of determining the total dust mass is insensitive to any cold dust — undetectable with 
Spitzer — but it is difficult to imagine how a significant fraction of interstellar dust (consisting 
almost entirely of small grains, a < 0.5 /xm) could remain cold within the hot plasma behind 
such a fast shock. Our value for Mjust includes only the dust now radiating along the NW 
filament; the overall faintness of the rest of the SNR at 24 /xm, combined with the lack of 
detection at 70 /xm, precludes determination of a dust mass for the entire remnant. 



5.3. The Ambient Dust-to-Gas Mass Ratio 



In previous works ( Borkowski et al. 2006 : Sankrit et al. 2010 : Williams et al. 2011bl lal) 
we have used the density determination from IR data, combined with the X-ray emission 
measure, to obtain a measure of the amount of shocked gas present within a given emit- 
ting volume. We have obtained the dust mass from the overall normalization to the dust 
spectrum, so dividing the two gives the dust-to-gas mass ratio. In the case of SN 1006, how- 
ever, the IRS spectra do not strongly constrain the density. Thus, we employ an alternative 
approach. 



Laming et al.l (119961 ) examined the ultraviolet spectrum of a section of the NW filament 
obtained with the Hopkins Ultraviolet Telescope (HUT) and measured 0.02 Hell (1640 A) 
photons cm~^ s~^. More recent work by Chayer (private communication) has revised the 
reddening towards SN1006 to be E{B — V)= 0.09, lowering the Hell fiux to 0.017 photons 
cm~^ s~^. According to the Laming et al. models, 0.0067 Hell photons come out for each H 
atom crossing the shock. For the 200" section of the filament observed (corresponding to a 
distance of 6.5 x 10^® cm at d = 2.2 kpc), the He II flux measured at Earth would be 



"Hell 



0.0067 (6.5 X 10- 

H atom J ^ 



cm 



And? 



where hq is the pre-shock H number density, Vg is the shock speed, and / is the line-of-sight 
depth through the emitting material, all of which are presumed (for this purpose) to be 



0.017 photons cm 



we can rearrange the above equation to 



unknown. Taking Fneii 
give UqVsI = 2.26 x 10^^ cm^^ s^^. The length of the NW-1 region is 3', or 5.9 x 10^ cm 
aX d = 2.2 kpc, so 1.3 x 10^^ H atoms per second enter the shock. If we assu me a standard 
post-s hock flow of 70 mas yr^^, based on the proper motion measurements of I Winkler et al. 



(120031 ) . then the width of the IR fllament corresponds to a timescale of 150 years. During 
this time, the shock will have swept up 6.3 xlO^^ H atoms, or 7.6 xlO^^M© of gas. The 
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uncertainties in the HUT measurement of the He II flux and in the reddening correction are 
about 15% each, while the uncertainty in proper motion is only 3% and the uncertainty in 
the relative atomic rates is probably about 30%, so we estimate an uncertainty in the gas 
mass of ~ 40%. Using the dust mass determined above, this leads to a dust-to-gas mass 
ratio of ~ (1 .5 ± 1.0) x 10~'^. The exp ected dust-to-gas mass ratio for the Galaxy within 
the model of iMazzei fc Barbard (120111 ) is 5.3 xlO~^, a factor of ~ 3 higher than what we 
measure. 

The dust-to-gas mass ratio in the ambient medium surrounding SN 1006 appears to be 
lower than that expected from dust models for the Galaxy, including those at high Galac- 
tic la titudes. This ha s been the case for virtually all SNR s studied, both in the Galaxy 



(e.g.. iBlair et al 



( iBorkowski et al. 



2007 : 



2006 



Arendt et al 



2010 



Williams et al. 



2006 



Lee et al.l |2009| ) and in the Magellanic Clouds 
2011bh . The ultimate reason for this remains 



a mystery, but in the case of SN 1006 specifically, its location within the Galaxy may pro- 
vide a clue. At some 500 pc above the Galactic p lane, SN 1006 is located at the interface 
between the disk and halo (see lPutman et al.ll2012l . for a review of gaseous galaxy halos and 
their interfaces to galactic disks ") . Accumulating evide nce points to a lower dust content 
there than in the Galactic disk (iBen Bekhti et al.ll2012l) . coupled with a higher proportion 
of smaller grains (jPlanck Collaboration et al.l 1201 if ). An expected increase in the sputtering 
rate for smaller grains might better match the rapid drop-off in the 24 /im profile behind the 
shock (Fig. [6]) than the current dust model, and it would also raise our i nferred value of the 
pre-sh ock dust-to-gas mass ratio, bringing it even more in line with the iMazzei fc Barbaro 
(I2OIII ) value. We further note that any dust that is too cold to be observed with Spitzer 
would increase the dust-to-gas mass ratio above our inferred value. 

Our current determination of the dust-to-gas mass ratio is quite uncertain, as we do not 
know the grain size distribution in the pre-shock gas. It is possible that models with a larger 
dust-to-gas ratio and a larger overabundance of small grains can be as viable as our current 
model with Ry = 2.0. This model degeneracy might be broken by future joint IR, X-ray, 
and UV studies focused on accurate determination of abundances of refractory elements in 
the postshock gas. 



6. Summary 

S'pzfeer observations of the entire region surrounding the remnant of SN 1006 show a clear 
detection of the NW filament at 24 fim, the same filament that is prominent in Ha emission. 
IRS spectroscopy of this filament shows that the IR is produced entirely by continuum 
emission from warm dust grains, with a spectrum peaking somewhere beyond 35 /xm. The 
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filament is not detected at 70 //m, likely due to the significantly higher background IR 
emission. The IR emission from the NW filament reaches its peak just inside of the Ha 
emission, while the X-ray emission peaks still a little further behind. This is as expected, 
since Ha emission should trace the exact location of the forward shock, while dust grains 
require a short amount of time to be heated to temperatures observable by Spitzer and then 
are sputtered away, while the post-shock gas requires even longer to become ionized enough 
to emit strongly at Chandra energies. 

Fits to the Spitzer IRS spectrum taken at two positions along the NW filament, using 
models of collisionally-heated dust grains, do not tightly constrain the post-shock density 
Hp. The nominal best fit gives Up = 1.4 ± 0.5 cm^^ (90%-confidence); we have adopted a 
model with rip — 1 cm~^, the value suggested from observations in other bands, for our 
subsequent analysis. The decay of the IR spatial profile behind the shock is far more rapid 
than expected from standard dust models. This is may be due in part to a high proportion 
of small dust grains in the ambient ISM at the location of SN 1006, high above the Galactic 
plane, and probably also results in part as a relic of the shock's having encountered a higher 
density region to the NW relatively recently. Nevertheless, it appears that sputtering rates 
significantly higher than those predicted by standard models will be required to fully explain 
the rapid decline in IR emission behind the shock. 

Calculations of the dust-to-gas mass ratio in the ambient ISM indicate a value of ap- 
proximately 1.5 X 10"'^, compared to the expected value of 5 xlO^'^ for the Galaxy at the 
latitude of SN 1006. This difference of a factor of several is a phenomenon seen in virtually 
all IR/X-ray analyses of SNRs. The remnant of SN 1006 provides a rare opportunity to study 
the ISM high above the Galactic plane, and further IR, X-ray, and UV studies are necessary 
to understand this remnant. 

This research is based on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology under a 
contract with NASA. We are grateful to the support staff at IPAC for their guidance with 
some of the subtleties of Spitzer data analysis. We also acknowledge thoughtful comments 
from the anonymous referee, which have prompted us to, we hope, clarify much of this 
paper. Primary financial support for this project has been provided by NASA through RSA 
1330031. PFW acknowledges additional support from the NSF through grant AST-0908566, 
and KJB acknowledges additional support from NASA through grant NNX11AB14G. 
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Fig. 1. — Images of SN 1006 in various bands. Panel a: Mosaic image at 24 /xm, assembled 
from all four MIPS AORs. b: The same 24 /xm mosaic, after subtraction of most point 
sources, displayed with a harder stretch, c: Deep Ha image, after continuum subtraction to 
reveal the extremely faint shell of emission to the S and E (the NW filament is saturated in 
this display; the irregular donut -like patterns are artifacts, the residual images from bright 
stars — from lWinkler et al.ll2003l ). d: M osaic image at 70 /xm, f rom a ll four MIPS AORs. e: 
Chandra AClS-l mosaic of SN 1006, from lCassam-Chenai' et al.l ( l2008l ): this image emphasizes 
the synchr o tron- dominated emission in the NE and SW. /.■ 1.4 GHz radio image, from 
Dyer et al.l (120091 1. All the images are aligned, with a field 40' square and oriented north up, 
east left. 



Fig. 2. — This 3-color image of the NW portion of SN 1006 shows the evolution in conditions 
behind the s hock, which is propa gating toward the NW (upper right corner). Ha filaments 
(green, from lWinkler et al.ll2003l ) form immediately behind the shock; 24 /im emission (red) 
from warm dust trails slightly b ehind as the dust grains are first heated and then destroyed; 
X-rays (blue, 0.5 - 2.0 keV, from llong et~aDl2003h peak further behind due to the ionization 
time to reach the species that are most effective at producing X-ray emission. The field is 
6.5' square; yellow regions indicate locations for the profiles shown in Fig. 3. 
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Fig. 3. — Radial profiles of emission in Ha, 24 /xm, and X-rays, taken perpendicular to the 
shock, for the four regions shown in Figure |2l (The data sources are the same as for that 
figure.) In each case, the profiles run from outside the shell inward, and the shock location 
is arbitrarily taken as the point at which Ha emission first rises sharply. 
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Fig. 4. — IRS positions superposed on the 24 /xm image of the NW filament of SN 1006. Red 
rectangles correspond to the second order (longer wavelength) spectrographs, while yellow 
corresponds to first order. The three numbered slit positions are from long pointings with 
the IRS Long-Low instrument (14-40 /xm). As we discuss in the text, we use positions 1 
and 2 for our spectral modeling. The shorter, nearly vertical slit position is from the IRS 
Short-Low instrument, where no emission was detected. 
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Fig. 5. — IRS spectrum, in black, combining positions 1 and 2 from Figure HJ and containing 
both orders of the Long-Low spectrograph. The points shown here, at wavelengths > 21/im, 
are all from order 1. The blue curve shows the nominal best-fit model, and gives a post- 
shock density in the filament of Up = lA cm~^. The red curve shows a fit with the density 
constrained at Up = 1.0 cm^^, within the limits of acceptable fits, and consistent with 
observations in other bands (text. Section 5.1). 
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Fig. 6. — A comparison of radial profiles of the 24 /xm emission (same as Fig. [3], scaled to a 
common peak intensity) with a model calculation, which has the post-shock proton density 
fixed at Up = 1 cm~^, consistent with observations in other bands (text, Section 5.1). A 
model assuming expansion into a uniform medium and normal sputtering rates decays far 
more slowly than the observations. 



